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Right-Handed Sneutrino as Cold Dark Matter of the Universe 
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We consider the minimal supersymmetric standard model (MSSM) extended by introducing three 
right-handed (s)neutrinos to account for neutrino masses in the oscillation experiments. Assuming 
that the neutrino masses are purely Dirac-type, the lightest right-handed sneutrino vr can be the 
lightest superparticle (LSP), which is a good candidate of cold dark matter (CDM) of the universe. 
We study the possibility of realizing j>_r-CDM, paying a special attention to the production of vr 
via decay of the next-to-lightest superparticle (NLSP) after its freeze-out time. It is shown that 
the late decay of the MSSM-LSP (the LSP among superparticles in the MSSM) can produce a 
sufficient amount of vr to explain the observed dark-matter density, and that the £r-CDM scenario 
can be realized in a wide range of parameter space. We also consider the constraint on the decay of 
MSSM-LSP from the big-bang nucleosynthesis (BBN); we found that the case with stau being the 
MSSM-LSP is severely constrained. 

PACS numbers: 14.60.Pq, 12.60.Jv, 95.35.+d, 98.80.Cq 



I. INTRODUCTION 

Supersymmetry is one of the most attractive candi- 
dates of physics beyond the standard model since it may 
solve various open questions in particle physics. It can 
solve the hierarchy and naturalness problems and also 
it can realize gauge-coupling unification. Supersymmet- 
ric theories have a great advantage in cosmology as well. 
This is because it can provide a viable candidate of cold 
dark matter (CDM) of the universe which cannot be un- 
derstood in the framework of the standard model of parti- 
cle physics. The lightest supersymmetric particle (LSP) 
becomes stable with i?-parity conservation and can be 
CDM. Recently, the relic density of CDM in the present 
universe has been precisely determined by the WMAP 
observation m: 



fl DM h 2 = 0.105 



+0.007 
0.013 ' 



(1) 



where h ~ 0.73 [l[ is the present Hubble constant in 
units of lOOkm/sec/Mpc. Any candidate of CDM must 
explain this dark matter density. Crucial questions left 
to us are, then, what is the LSP for CDM and how they 
are produced accounting for ((T|) in the history of the uni- 
verse. So far, various scenarios have been discussed in 
literatures 0] . 

In Ref. [3j , we have proposed a scenario that the light- 
est right-handed sneutrino vr is the LSP and is CDM 
of the universe. 1 This is motivated by very small but 
non- vanishing neutrino masses strongly suggested by the 
experimental evidences of neutrino oscillations (see, e.g., 
Refs. d, 0, H, H, EH)- The non-zero neutrino masses re- 
quire physics beyond the standard model. The simplest 
way to generate neutrino masses is probably to introduce 



1 The possibility of the right-handed sneutrino as CDM has also 
been discussed recently in different context. See Refs. 



right-handed neutrinos (as well as sneutrinos in super- 
symmetric theories). 

It has been widely discussed that right-handed 
(s)neutrinos are introduced together with their super- 
heavy Majorana masses which are much larger than the 
electroweak scale ~ 100 GeV. Then the smallness of neu- 
trino masses is naturally explained by the so-called see- 
saw mechanism 11]. In this case, the right-handed sneu- 
trino cannot be the LSP. 

On the other hand, neutrinos can obtain very small 
masses without invoking the seesaw mechanism. By in- 
troducing right-handed (s)neutrinos with vanishing Ma- 
jorana masses, neutrinos become massive Dirac fermions 
after the electroweak symmetry breaking. In this case, 
as we will show later, the neutrino Yukawa coupling con- 
stants are of the order of 0(1O -13 — 10~ 12 ) or smaller in 
order to account for the neutrino mass scales in the os- 
cillation experiments. It should be noted that such small 
coupling constants are natural in the 'tHooft's sense 
This is owing to the fact that the chiral symmetry in neu- 
trino sector is recovered in the limit of vanishing neutrino 
Yukawa coupling constants. 

In supersymmetric theories, when neutrino masses 
are purely Dirac-type, right-handed sneutrinos receive 
masses dominantly from the effects of supersymme- 
try breaking. Throughout this paper, we consider 
gravity-mediation type models of supersymmetry break- 
ing. Then, all the supersymmetry breaking masses are of 
the order of 0.1-1 TeV. As a result, there is a possibility 
that the lightest right-handed sneutrino vr is the LSP. In 
Ref. 3, we have shown that the LSP vr is a good can- 
didate of CDM because it is stable and charge-neutral, 
and also because its relic density can become consistent 
with the dark matter density with relevant choice of pa- 
rameters. 

In the early universe, vr is not thermalized since its 
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interaction is extremely weak. 2 Thus, assuming that the 
initial abundance of vr is zero, 3 Vr should be produced in 
some processes after inflation if vr-CDM is realized. One 
possible production is via decays of superparticles [||. 
Superparticles in the minimal supersymmetric standard 
model (MSSM) are in chemical equilibrium when the tem- 
perature is high enough, so they are quite abundant in 
the early universe. Decays of those MSSM superparticles 
can produce vr. The decay processes may occur when 
the MSSM superparticles are in the chemical equilibrium 
and when they freeze-out from the thermal bath. Impor- 
tantly, the decay processes after the freeze-out may be 
important for the production of vr. This happens for 
the decay of the lightest superparticle in the MSSM sec- 
tor, which we call the MSSM-LSP. This is because the 
neutrino Yukawa coupling constants are so small that 
the lifetime of the NLSP becomes rather long. Such late 
decay of the NLSP into vr can be one important source 
of ^-CDM. 

In this paper, we consider the scenario with C^-LSP 
and investigate its production in the early universe. In 
Ref. [3|, we have mainly analyzed the production by de- 
cays of superparticles which are in chemical equilibrium. 
Here, we study another production process in detail, that 
is, decay of the MSSM-LSP after its freeze-out time. We 
estimate the relic density of vr from such a decay adopt- 
ing the minimal supergravity model (l3j . and also dis- 
cuss the implication of the vr-CDM scenario. We also 
reconsider the decay of MSSM superparticles in chemical 
equilibrium. 

The outline of this article is the following. We start 
by reviewing the framework of our analysis in Sec. [Til 
where we present properties of neutrinos and sneutrinos 
under consideration. In Sec. IIIII various production pro- 
cesses of vr are discussed. We hrst study the production 
of vr by decays of the NLSP after the freeze-out time. 
The constraint on the NLSP decay coming from the big- 
bang nucleosynthesis (BBN) is also discussed. We then 
turn to reconsider the production by decays of superpar- 
ticle in chemical equilibrium to complete the discussion 
in Ref. In Sec. IIIII other possible mechanisms of 

the production of vr are also mentioned. Finally, we 
conclude in Sec. IIVI 



II. FRAMEWORK 

First of all, let us explain the framework of our anal- 
ysis. We consider the MSSM with three generations of 
right-handed (s)neutrinos, where neutrino masses are as- 
sumed to be purely Dirac-type. The superpotential is 



then given by 

W = y v H u -Lv R -y e H d -L t R + n H H d -H u , (2) 

where we have omitted terms with quark superficlds 
since they are irrelevant for our discussion. Here H u — 
(H+,H°) and H d = (H%,HJ) are the Higgs superfields 

coupled to up- and down-type quarks, and L = {v^,^) 
are the left-handed lepton superfields. The superfields 
of right-handed neutrinos are denoted as vr. (In this 
article, "hat" is for superfields while "tilde" is for super- 
particles with odd i?-parity.) fin is the supersymmetry- 
invariant Higgs mass. Here and hereafter, the generation 
indices are implicit for simplicity. 

In this framework, neutrinos obtain masses as 

m v = y v (H°) = y v v sin/3 , (3) 

where v ~ 174 GeV is the vacuum expectation value 
(VEV) of the standard-model-like Higgs boson, and 
tan/3 = Numerically, we obtain 

y v sin/3. 3.0 x 10- x . (4) 

Importantly, neutrino-oscillation experiments have pro- 
vided only the mass-squared differences of neutrinos [1, 

Nat„^ 2 - 8Xr3eV2 . ( 5 ) 

K] s „, ar ^.9xl0^eV 2 , (6) 

and the absolute scales of neutrino masses have not been 
determined. There are then three possible mass spectra: 

(i) the normal hierarchy case, (ii) the inverted hierarchy 
case, and (hi) the degenerate case. On the other hand, 
the cosmological observations place the upper bound on 
m„. Here, we adopt the upper bound derived from the 
WMAP three-year data: Y, m » < 2 -° eV ( 95% CL ) R 
where summation is over all the three neutrino flavors. 4 
Then, we obtain 

m v < 0.67 eV. (7) 

In the case with hierarchical neutrino masses (i) and 

(ii) , the heaviest neutrino mass is well approximated 
by ^/[Am 2 ] atom and the corresponding neutrino Yukawa 
coupling constant, which is the largest one, is about 
j/ySin/3 ~ 3.0 x 10~ 13 . On the other hand, in the de- 
generate case (hi), the coupling constants become larger; 
if the neutrino masses take their largest possible value 



2 If it was thermalized by some unknown interactions in the very 
beginnings of the universe, the present relic density of vr would 
overdose the universe and be inconsistent with the observation. 

3 See, however, the discussion in Sec. IIII Dl 



The inclusion of the data other than the cosmic microwave 
background radiation (CMBR) observation (e.g., the large scale 
structure and the matter power spectrum inferred from the Ly- 
ct) makes the mass limit of neutrinos more stringent (see, for 
example, Ref. p]]). We take here a conservative approach and 
use the bound only from the CMBR data. 
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of ~ 0.67 eV, Yukawa coupling constant is as large as 
y v sin/3 ~ 3.8 x 10~ 12 . We can see that the neutrino 
Yukawa coupling constants are very small in any case, 
which are much smaller than those of other quarks and 
lcptons. 

Next, we turn to discuss sneutrinos. For this purpose, 
let us introduce the soft supersymmetry breaking terms 
with sneutrinos as 

-£ soft D M 2 \L\ 2 + M 2 Ji> R \ 2 

+ (X H u -LD R -A e H d -L£ e R + /i.e.) ,(8) 

where Ml, M Vr and A v are supersymmetry breaking 
mass parameters. We parameterize A v as 

A v =y v A v = y v a v M L , (9) 

where a v is a dimension- less constant and |et„| is expected 
to be of the order of unity in simple models of super- 
gravity. Due to the smallness of the neutrino Yukawa 
coupling constants, the left-right mixing of sneutrinos is 
very small. Thus, we can treat the left- and right-handed 
sneutrinos as mass eigenstates. Their masses are given 
by 

mj L -Ml + \ cos 2/3 m| , m? fi ~ M 2 R , (10) 

where we have omitted the negligible contributions of 
m 2 . The left-right mixing angle of sneutrinos is denoted 
by 9, which is found as 

tan26= 2m - |c ° t ^V^ 1 , (11) 

and O is highly suppressed because of the smallness of 
m v . It can be seen from (fit)]) that masses of right-handed 
sneutrinos are determined solely by the supersymmetry 
breaking masses, and hence the lightest right-handed 
sneutrino may become the LSP. 

III. PRODUCTION OF THE LSP 
RIGHT-HANDED SNEUTRINO 

As we have shown in Ref. [3], the LSP right-handed 
sneutrino is a viable candidate of CDM. This is be- 
cause it is stable under the i?-parity conservation, and 
also because it has only suppressed interactions, i.e., the 
neutrino Yukawa interactions and the gauge interactions 
through the left-right mixing. 

In this section, we discuss production processes of 
right-handed sneutrino v R in the early universe and cal- 
culate the density parameter of right-handed sneutrino: 

n, R = ^, (12) 

Per 

where Pv R ,<a is the present energy density of v R and 
p cr = 1.05 x 10 -5 h 2 GeV/cm 3 is the critical density. Im- 
portantly, there are two types of production processes of 



right-handed-sneutrino LSP; one is the decay of MSSM 
superparticles in the chemical equilibrium and the other 
is the decay of the MSSM-LSP after freeze out. We have 
considered the former effect in Ref. Q and found that, in 
most of the cases, Q,p R via the decay of MSSM superpar- 
ticles in the chemical equilibrium becomes smaller than 
Hdm- Thus, it is probable that the relic v R are mostly 
from the decay of the MSSM-LSP after freeze-out. So, 
we first consider the latter effect in detail although we 
will also discuss the former production process later. 



A. D R production via MSSM-LSP decay 

Let us discuss the production of v R via decays of the 
MSSM-LSP after freeze out. As we will discuss, the den- 
sity parameter of v R is sensitive to the MSSM parame- 
ters. In this paper, we will not study the complete param- 
eter space since it is beyond the scope of our purpose; we 
rather illustrate the basic features using a simple model 
of supersymmetry breaking. Here we consider, as an ex- 
ample, the minimal supergravity model with the gauge 
coupling unification at Ma — 2 x 10 16 GeV. This model is 
described by mo (the universal scalar mass at Mq), m 1( / 2 
(the unified gaugino mass at Mq), and Aq (the universal 
A-parameter at Mq) together with tan/3 and the sign of 
fin ■ Then, the on-shcll mass of v R is given by 

mo R = mo , (13) 

since the renormalization group (RG) evolution of m^ R 
is negligible due to the smallness of y v . 

Adopting the minimal supergravity model, three right- 
handed sneutrinos are almost degenerate. The mass dif- 
ferences are mainly from the effects of RG evolution 
through y v , and hence they are very small. 5 Conse- 
quently, all the three right-handed sneutrinos are stable 
within the age of the universe and contribute to the dark- 
matter density. 

Since the neutrino Yukawa coupling constants are neg- 
ligibly small, the decay of the MSSM-LSP into v R is ex- 
tremely suppressed. Detailed value of the lifetime of the 
MSSM-LSP depends on MSSM parameters. However, 
as we will show later, lifetime of the MSSM-LSP is long 
enough so that the decay of the MSSM-LSP occurs suf- 
ficiently after its freeze out. 

The density parameter of v R from the decay of the 
MSSM-LSP after freeze out is given by 

1l v r — "mssm-lsp , v 14 v 

TOmssm-lsp 

where fi MSSM _ LSP is the (would-be) present density param- 
eter of the MSSM-LSP for the case where it is stable. 



5 Other sources of the mass differences are the neutrino masses 
and the contributions from the left-right mixing of sneutrinos. 
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FIG. 1: The parameter range to realize i/r-CDM in the minimal supergravity model with Ac = (shown by green "+" marks). 
The solid line represents mc, R = wimssm-lsp; vr becomes the LSP below this line. The MSSM-LSP is the lightest neutralino 
X? or the lighter stau f\ in the region above or below the dashed line. The vertical dot-dashed lines show the lower bounds on 
m i/2 from the Higgs mass and b — > s*y. 



^mssm-lsp is estimated using conventional method since 
the neutrino Yukawa coupling constants are very small. 6 
vr can account for the present dark-matter density if 
fi|° = fi DM . We note here that flp° is insensitive to the 
reheating temperature Tr after inflation as long as Tr is 
higher than Tp, where Tp is the freeze-out temperature 
which is roughly given by Tp ~ Wmssm lsp /20. 

Based on Eq. (fT4"|) , z/r-CDM is realized with some 
choice of MSSM parameters which satisfies the following 
relation: 

WmSSM-LSP q -\ 
^mssm-lsp — "dm ■ k^*- 1 ) 



In this analysis, we assume that only the thermal relic of the 
MSSM-LSP contributes to Hmssm-lsp- The MSSM-LSP might 
be also produced via non-thermal processes [l5| . Even in this 
case, if Hmssm-lsp is given as H15I I. the fjj-CDM is realized. 



Since m M ssMLSP 7^ m v R i this implies that the MSSM 
parameters realizing j>r-CDM are different from those 
for the conventional scenario where the MSSM-LSP, say 
the lightest neutralino, becomes CDM. Furthermore, vr- 
CDM is possible even if the MSSM-LSP is an electrically 
charged and/or colored superparticle, and hence the vr- 
CDM is realized in a wider parameter range compared 
with the MSSM-LSP dark matter. 

In the parameter range of interest, the MSSM-LSP 
is the lightest neutralino x1 which is Bino-like or the 
lighter stau fi. We estimate their (would-be) relic den- 
sity r2 MSSM _ LSP /i 2 by using the micrOMEGAs package 
and then calculate ~° h 2 using Eq. (fl"4f . In Fig. [TJ we 
show the parameter region where fi|°/i 2 is in the range 
((T|). Here, for simplicity, we set Aq = 0. In the figure, we 
also impose the following phenomenological constraints; 
the mass of the lighter neutral Higgs boson (the standard- 
model like Higgs boson) be larger than the experimental 
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FIG. 2: The mass of i/R-CDM in the minimal supergravity 
model. We take here Aq ~ 0, tan (3 = 10 and \ih > 0. The 
thick line corresponds to Q^h 2 = 0.105, and the thin lines 
correspond to f^°/i 2 in the range (TTJ). The solid or dashed 
lines are for the Xi or f\ MSSM-LSP, respectively. 



bound llj 



m h o > 114.4 GeV, 



(16) 



much larger than ~ 1 TeV. For m^o 3> 700 GeV, the 
yield variable of fi , which is defined \>y the ratio of the 
number density to the entropy density, is given by [l9[ 



10" 12 X Cy 



1 TeV 



(18) 



where cy is a constant of the order of unity; we numer- 
ically found that cy = 1.2 — 1.3 in the parameter region 
we are interested in. We can then obtain the relic density 
of vr from fi decay as 



0.027 x cy 



VI TeV/ V 



100 GeV 



(19) 



The numerical calculation shows that, when mo R = mo = 
100 GeV, for example, ^l° R h 2 = 0.105 is realized for 
rrifx — 3.1 TeV, which corresponds to m 1 / 2 — 8.7 TeV 
and m^o ~ 4.1 TeV. Even when we take larger value 
of mfj (and m^o), the ur-CDM is possible by taking 
m v R ^5 100 GeV, as seen from Eq. (fl9|) . 

Now, we come to a position to discuss the lifetime of 
the MSSM-LSP. As we mentioned, Eq. (H} is applicable 
when the MSSM-LSPs decay after the time of freeze out. 
In order to see if this is the case, let us estimate the de- 
cay rate of the MSSM-LSP. When the MSSM-LSP is the 
Bino-like neutralino, it decays into vr and anti-neutrino 
(and also its CP conjugate state) through the left-right 
mixing of sneutrinos. The decay rate is estimated as 



and the branching ratio of b — > sj be in the three-sigma 
range of the observational data [18| : 



Br(& -> sj) = (3.55 ± 0.78) 



x 10 



-4 



(17) 



One can see that vr can be CDM of the universe in a 
wide parameter region. In addition, we emphasize that 
vr-GDM can be realized not only in the case where \i 1S 
the MSSM-LSP but also in the case with ^-MSSM-LSP. 
This is an interesting feature of the present scenario. We 
also stress here that is insensitive to the neutrino 
Yukawa coupling constants. 

In Fig. [21 we show the mass of right-handed sneutrino 
which gives the correct value of f2~° to realize ^-CDM 
as a function of the mass of the MSSM-LSP. Here, we 
take Aq = 0, tan j3 = 1Q and hh > 0; in this case, when 
n F ~°h 2 = 0.105, the MSSM-LSP is the lightest neutralino 
for m M ssM-LSP < 510 GeV, and is stau for to M ssm-lsp ^ 
510 GeV. Fig.[2]shows that the mass of the right-handed 
sneutrino is bounded from above to realize ^-CDM. 

In the region where Xi becomes the MSSM-LSP, we 
find an approximate relation m^o ~ 1.8mj H to realize 
= Q, DU . In this case, the mass of the 9r-CDM is 
well determined from rn^o even if one varies f2?° h 2 in the 
range given in Eq. (fTJ) ■ Moreover, we find that masses of 
Xi and ti are close to each other and their co-annihilation 
effects is important in estimating f2 MSSM _ LSP . (See the 
following discussion and Fig.|4j) 

As m^o becomes larger, the MSSM-LSP is f\. Then, 
the relation f2?° = f2 DM is realized even when m^.o is 



r 



"32^ m *? 1 



2 



(20) 



where g\ is a gauge coupling constant of U(l)y. No- 
tice that the Bino-like neutralino Xi universally couples 
to (s)neutrinos in all three generations. Thus, the de- 
cay process into the sneutrino which interacts via the 
strongest Yukawa interaction is the most important since 
the mixing angle O is proportional to the Yukawa cou- 
pling constant. 

When the MSSM-LSP is n, the situation is slightly 
more complicated since the dominant decay mode de- 
pends on the mass difference between ?i and v. When 
mfi > TUu B + mw , the decay rate is given by 



>v R W 



gl Q 2 itt(ti) |2 |ttMNS|2 m n 



1 - 



2(m? 



2 

m% 



+ 



w 



i 



mi 



3/2 



,(21) 



where g 2 is a gauge coupling constant of SU(2)j,, m\y 
is the H^-boson mass, and is the mixing matrix 

of staus which relates the gauge eigenstates and mass 
eigenstates as 
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(22) 
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FIG. 3: The lifetime of the MSSM-LSP yielding n~°/i 2 = 
0.105 in the minimal supergravity model. We take here Ac ~ 
0, tan/3 = 10 and fi H > 0. The MSSM-LSP is Xi or fi 
when m-MSSM-LSP is smaller or larger than about 510 GeV. 
The solid lines are the lifetime assuming the normal hierarchy 
of neutrino masses and taking m v — J[Am'] Boa , while the 
dashed lines are the lifetime assuming the degenerate neutrino 
masses and taking m„ = 0.67 eV. 



where f\ and f 2 are mass eigenstates with < m^ 2 . 
In addition, [/ MNS is the neutrino mixing matrix. Here, 
we use 

jjMNS = 1 /^/2_ Notice that the aboye 

expres- 
sion is relevant for the case with hierarchical neutrino 
mass; for the degenerate case, decay processes into all 
three sneutrinos may be equally important. For exam- 
ple, when the mass differences among three left-handed 
sleptons are negligible, |[/^ NS | 2 should be replaced by 
1. Note that f\ may also decay into vr and the charged 
Higgs boson. Such a decay channel is, however, kine- 
matically forbidden in the parameter region which we 
are interested in. When the decay process fi — * vrW is 
kinematically forbidden, t\ mainly decays as f\ — > Ur£v, 
vr q q' with exchanging virtual M^-boson. 

Using y v ~ 3.0 x 10~ 13 , which is the suggested value of 
the neutrino Yukawa coupling constant in the case with 
the normal hierarchy neutrino mass matrix, the lifetime 
of the MSSM-LSP for the case yielding ^h 2 = 0.105 
is shown in Fig. [3l Here, we take Aq — to illustrate 
the typical value of the lifetime. We can see that the 
lifetime is roughly of the order of 10 2 -10 3 sec, and hence 
the decay of the MSSM-LSP occurs after its freeze-out 
time. When neutrino masses are degenerate, the lifetime 
of the MSSM-LSP becomes significantly shorter because 
the neutrino Yukawa coupling constants are enhanced 
and also because the MSSM-LSP decays into three (al- 
most) degenerate right-handed sneutrinos. Even in this 
case, the MSSM-LSP decays into vr after the time of 
the freeze-out. Therefore, calculation with Eq. ([H)) is 
justified with any of the neutrino mass matrix. 

Here, we should comment on the fact that the lifetime 



of the MSSM-LSP depends on the A-parameter since the 
left-right mixing angle of the neutrino is sensitive to it. 
By solving RG equations, we found that A v parameter 
at the scale mz is given by 



A„ ~ A ( 



0.59 m 1 / 2 



(23) 



Thus, as A G changes, the lifetime of the MSSM-LSP 
varies. 

We would like to note that the vr-CDM scenario and 
MSSM-LSP dark matter scenario are realized in differ- 
ent parameter space. In particular, we emphasize that 

^MSSM-LSP 

M is required if the right-handed sneu- 
trino produced from the decay of MSSM-LSP after freeze 
out becomes CDM. Using this fact, we can potentially 
distinguish the scenario of Dr-CBM and MSSM-LSP 
dark matter with precise determinations of the properties 
of MSSM particles by future collider experiments. To see 
this, in Fig.|H we show the regions on to^o vs. —m^o 
plane where the density parameter of the LSP becomes 
consistent with the presently observed dark matter den- 
sity for the cases where the LSP is vr and \v We can 
see that, for a given value of m^o, the suggested values 
of are different by 5 GeV or so between two scenar- 
ios. We also show the correlations between mf, and m~+ 

Xi 

(with X\ being the lighter chargino) in two scenarios. 
We can see that, for a given m^, the suggested value of 
m-+ is different by ~ 10 GeV. Thus, with the precise 
measurements of the masses of MSSM superparticles by 
future collider experiments, like the LHC and ILC, we 
may be able to distinguish the vjj-CDM and the MSSM- 
LSP-CDM scenarios. 



B. Constraints from BBN 

As one can see from Fig. |3l the lifetime of the MSSM- 
LSP may be so long that the decay of the MSSM-LSP 
occurs around or after the epoch of BBN. If so, abun- 
dances of light elements which are produced by the stan- 
dard BBN reactions are affected by the decay of the 
MSSM-LSP and, consequently, the success of BBN may 
be spoiled. 

Effects of late-decaying particles on BBN are inten- 
sively studied in Ref. [20]; the most important effects 
of the late-decaying particles are from p <-* n conver- 
sion, hadro-dissociation, and/or photo-dissociation pro- 
cess, depending on the lifetime. For the case of our in- 
terest, lifetime of the MSSM-LSP is < 10 4 sec. In such 
a case, p <-> n conversion and hadro-dissociation are im- 
portant. If the abundance of the MSSM-LSP is too large, 
the abundances of light elements are too much affected 
to be consistent with observations by these processes. In 
order not to spoil the success of BBN, upper bound on 
the combination Bhad-E'vis^ssM-LSP is obtained. In our 
case, .Bhad is the hadronic branching ratio of the MSSM- 
LSP, E v i s is the net energy carried away by hadrons in 
the decay of the MSSM-LSP, and Y MSSM _ LSP is the yield 
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FIG. 4: The mass differences m?, — m-o and 2m f, — m~+ 
realizing £/r-CDM in the minimal supergravity model. We 
take here Ac = 0, tan (5 = 10 and fiH > 0. The thick solid 
line corresponds to Q^°h 2 = 0.105 and the thin solid lines 
correspond to f2~°/i 2 in the range JT]). For comparison, we 
also show the mass differences in the case where the MSSM- 
LSP Xl becomes dark matter by dashed lines. 

of the MSSM-LSP. In order to see if %-CDM is viable, 
we will calculate these quantities in the following. 

We first consider the case where the MSSM-LSP is 
Xi- Even though the dominant decay mode of Xi is 
into neutrino and right-handed sneutrino, there still exist 
hadronic decay modes: 

Xi^VRvqq, V* K vqq, v R i + qq , v* R f qq , (24) 

where q and q' are quarks and £ are charged leptons. 
These processes are mediated by on-shell and/or off-shell 
gauge bosons. We have numerically calculated the decay 
rate of these processes and obtained Shad and E v i s . Re- 
sults for the case with Aq — 0, tan/3 = 10 and \ih > 
are shown in Fig. As one can see, hadronic branching 
ratio becomes O(10~ 4 — 10 -3 ) when the mass of xi is rel- 
atively large. This is because, when m^o is large enough, 



FIG. 5: The hadronic branching ratio Bhad and the net energy 
carried by hadrons _E v is in the x° MSSM-LSP decay. Here 
we consider the minimal supergravity model with Ac = 0, 
tan/3 = 10 and [in > 0, and require Qj®h 2 — 0.105. 



the decay modes Xi ~^ v R l + W~ and Xi —> v rvZ (and 
their CP conjugated processes) are kinematically allowed 
and hence hadrons can be produced via the "three-body" 
decay processes. On the contrary, when m^o is small, the 
weak bosons are always off-shell and -Bhad is suppressed. 

For MSSM parameters which give correct value of 
0~° to realize i/jj-CDM, we calculate the product 
■Bhad-Evis^MssM-LSP as a function of the MSSM-LSP mass. 
The results are shown in Fig. [6] For the same parame- 
ter set, we also calculate the lifetime of the MSSM-LSP. 
Once the lifetime of the MSSM-LSP is given, we can es- 
timate upper bound on the product Bhad-Svis^ssM-L 
order not spoil the success of BBN. When the lifetime is 
longer than ~ 100 sec, the most important bounds come 
from the abundances of D/H and 6 Li/H; here we adopt 
constraints, which are conservative ones, based on Fig. 39 
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FIG. 6: The product Bhad-E'visi^ssM-LSP in terms of the 
MSSM-LSP mass in the minimal supergravity model, (shown 
by the solid lines). The upper or lower panel is for the case 
when Ac — or Ac = — 3mo, respectively. Here we take 
tan/3 = 10 and fi H > 0, and require Sll°h 2 = 0.105. The 
MSSM-LSP is Xi or Ti f° r w-mssm-lsp is greater or smaller 
than about 510 GeV. The dot-dashed lines show the upper 
bound from the BBN constraint on the MSSM-LSP decay by 
assuming the normal hierarchy of neutrino masses and taking 
m„ = ^/[AmJ] at om- The dashed lines show the bound by as- 
suming the degenerate neutrino masses and taking m„ = 0.67 
eV. 



in Ref. \2£ 




for 12 > x > 3.1 
for 3.1 > x > 2.6 
for 2.6 > x > 1.6 



(25) 



where x and y are 

X = l°glo( T MSSM-LSp/ SeC ) 5 



V = log'lO^had-Evis^ 



MSSM-LSP 



/ GeV) . 



(26) 



With shorter lifetime, overproduction of 4 He gives an up- 
per bound on the abundance of the MSSM-LSP; here we 



show the constraint based on the observational constraint 
given in Ref. [2lj : 



-0.17 a; -9.6 
y < \ -0.88x-9.4 
-3.8 a;- 10.6 



for 3.5 > x > 0.39 
for 0.39 > x > -0.41 
for -0.41 > x > -1.2 



(27) 



Upper bound on Bhad-Evis^MssM-LSP is also shown in 
Fig. [6] Notice that the lifetime of MSSM-LSP depends on 
the MSSM parameters and also on the neutrino Yukawa 
coupling constants, so we show the upper bound for sev- 
eral cases. 

With hierarchical neutrino mass matrix, we can see 
that i/R-CDM is viable in most of the cases, although 
some parameter region with large m^o may be con- 
strained from the BBN, in particular for the case with 
Aq = 0. If we consider the degenerate case with 
m„ = 0.67 eV, lifetime of Xx becomes much shorter and 
hence the bound becomes weaker. Therefore, f>^-CDM 
does not conflict with the BBN constraints irrespective 
of the value of Ac, as long as there is no cancellation in 
the low-energy value of A v (see Eq. ([2"3"])). 

Next, we consider the case where the MSSM-LSP is 
T\, As mentioned before, when fi is the MSSM-LSP, 
z/r-CDM is realized when mf x is relatively large. In this 
case, f\ mainly decays into vr and W. Since IF-boson 
decays into quarks with the branching ratio of 2/3, we 
estimate the hadronic branching ratio of f\ as 



B 



had 



In addition, £" v is is given by 



2m.f 



(28) 



(29) 



In Fig. El we also plot -Bhad-Evis ^mssm-lsp as a function 
of the mass of fi . Compared to the case with ^"MSSM- 
LSP, -Bhad-EvisYMssM-LSP becomes much larger since the 
hadronic branching ratio is close to 1 in this case. In par- 
ticular, with the normal hierarchy of the neutrino masses, 
lifetime of fi is of the order of 10 3 sec or longer as far 
as \a v \ ~ 0(1). In this case, the whole region of the 
fi MSSM-LSP yielding = fl DM is excluded. No- 
tice that, although we have shown only the results when 
m fl < 1 TeV, the n MSSM-LSP with m fl > 1 TeV is 
also excluded since the lifetime and E v [ s become larger. 

On the other hand, for the case of degenerate neutrino 
masses, there appear the regions in which the fi MSSM- 
LSP is consistent with the BBN. In such regions, the 
lifetime of fx is sufficiently short and it only receives the 
weaker constraint (|27[) . We find that the possible masses 
of fx are relatively small and > 1 TeV is excluded 
by the BBN even with degenerate neutrinos. 

We also find that the light fx is also excluded by the 
BBN when tan j3 is large. This is because, in such a case, 
the mass difference between Vr and f\ is smaller than 
mwj and the decay channel fi — > vr-\-W is kinematically 
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blocked. If this is the case, the lifetime becomes longer 
and more stringent constraints apply. 

From Fig. [6j we can learn that BBN constraint puts 
the upper bound on the mass of the MSSM-LSP. Interest- 
ingly, this bound excludes the possibility of the MSSM- 
LSP being very heavy if vr is CDM, say being heavier 
than of the order of TeV. 

Before closing this subsection, we should comment that 
some of the results obtained in the previous and this 
subsections are from the fact that we use the minimal 
supergravity model. In particular, in the minimal su- 
pergravity model, masses of right-handed sneutrinos are 
equal to mo- In this case, we obtained upper bound on 
the mass of the MSSM-LSP if v R is CDM. This result 
may not hold if we consider other types of models of su- 
persymmetry breaking. In particular, if we adopt smaller 
sneutrino masses, we can push up the mass of the MSSM- 
LSP. Furthermore, if the mass of the left-handed sneu- 
trino becomes close to m£> R , the left-right mixing angle 
of sneutrinos are enhanced. (See Eq. (TTTj) .) In this case, 
the lifetime of the MSSM-LSP becomes shorter and, as 
a results, the BBN constraint becomes weaker. In such a 
situation, the fi MSSM-LSP is possible in a wider param- 
eter space. Finally, we should mention that the MSSM- 
LSP can be a superparticle other than the Bino-like neu- 
tralino or the lighter stau. To realize the vr CDM, the 
MSSM-LSP may be electrically charged or colored su- 
perparticle, if its relic abundance satisfies (fl5|) and if its 
decay is consistent with the BBN. 



C. Decays of superparticles in chemical equilibrium 

In the previous subsection, we have assumed that the 
vr production from the decay of MSSM superparticles in 
chemical equilibrium is negligible. Indeed, this is the case 
in most of the parameter region. In some case, however, 
such a production mechanism also gives sizable amount 
of right-handed sneutrino. (We denote its contribution to 
the relic density by fiSf.) In this subsection, we consider 
f*? E . 

Vr 

f2p E is estimated by solving the Boltzmann equation 



dnp R 
dt 



3Hn,y R 



Cd 



ccay , 



(30) 



where 



is the number density of Vr. Denoting the 



distribution function of particle x in the chemical equi- 
librium as f x , the decay term is given by 



Cdccay ^ ^ 



d 3 k 



(27T) 



(31) 



where ^ x — m x j ' \Jk\ + ra x is the Lorentz factor, (25^ + 1) 
is the spin multiplicity of x, and (1± is the averaged 
final-state multiplicity factor for fixed value of initial- 
state momentum (with the positive and negative signs 



being for bosons and fermions, respectively.) In addi- 
tion, the summation is over all the possible production 
processes of right-handed sneutrino. Then, the relic den- 
sity is found as 



j^CE _ m VR n vn.,0 



(32) 



where np^o is the present value of i%o R - 

Solution to the above Boltzmann equation has been 
already studied in Ref. [|[; in Ref. [H, we found that 
f2~^ becomes much smaller than ri DM with neutrino mass 
matrix with normal hierarchy if tri-linear scalar coupling 
is small. Indeed, if A v = 0, the relevant processes of the 
vr production are the decay of Higgsinos, H — > vrT>l 
and H + — > vrI\. Then, we found 



< 1-7 x 10' 



2.8 x 10" 3 eW 



(33) 



(Here, we assume that the neutral and charged Higgsinos 
are mass eigenstates with mass |/iff|.) The above expres- 
sion is independent of m^ R , and the maximal abundance 
is obtained when ~ 2.75m^ R . 

It is notable that f2~™ can be as large as Hdm when 
the mass of left-handed sneutrino is close to right-handed 
sneutrino mass. This is because the left-right mixing 
angle of sneutrinos is enhanced in such a case. In the 
case of minimal supergravity model, such a mass spec- 
trum is hardly realized. In general, however, left-handed 
and right-handed sneutrino masses are free parameters 
and hence those masses may be degenerate. Thus, in 
this subsection, we work in the general framework of 
the MSSM rather than adopting minimal supergravity 
model. When the mass difference between left-handed 
and right-handed sneutrinos are small, the right-handed 
sneutrino production is dominated by the decay of Wino 
and Bino in the chemical equilibrium: W° — > vrVl, 
W + — > vrI\, and B — > vrv^ (and their CP conjugated 
processes). In the following, let us discuss this case ap- 
proximating that the mass eigenstates of charginos and 
neutralinos are Wino and Bino (and Higgsinos). 

The left-right mixing angle of sneutrinos in the vacuum 
is given in Eq. (jlip . Importantly, at high temperature, 
the mixing angle varies since the expectation value of the 
Higgs boson depends on the temperature; we estimate the 
temperature-dependent mixing angle as 



6(T) 



1 tan _i pi/.-rlDlcot. //«/, 



A* 



m% L (T) - ml r 



(34) 



where v(T) is the temperature-dependent expectation 
value of standard-model-like Higgs boson Hsm- Here, 
we assume that the heavier Higgses have masses much 
larger than mz and that the temperature dependence of 
Higgs mixing angle (3 is negligible. In addition, m~ L (T) 
is the temperature-dependent mass of left-handed sneu- 
trino. With 0(T), decay rates of the relevant processes 
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are given by 



^W°-> VRVL 



W+^vrIJ 



J3—*v R v L 



J-6 2 (T)m # [ 1 - 
64tt v ' w \ to 2 - 



2 

TO- 

VR 



mi 



w . 



64tt y ' b 



(35) 
(36) 
(37) 



In our analysis, we approximate the Higgs potential in 
thermal bath as 1221 



2\2 



+ — (2m 2 w + m 2 z + 2 m 2 )T 2 \H SM \ 2 . (38) 

Minimizing this potential, we obtain v(T); denoting crit- 
ical temperature (which is defined as the temperature 
where the curvature of Vr at Hsm — vanishes) as Tc, 
we obtain the expectation value at T < Tc as 



v(T) = vJi-Tyr c 



(39) 



175 



In our calculation, we use m,h = 115 GeV and mt 
GeV. In this case, T c = 139 GeV. 

With the shift of the expectation value of Higgs bo- 
son at high temperature, the mixing angle is modified, 
as shown in Eq. (|34|) . The effects are from the change 
of the left-right mixing mass of sneutrino, and also from 
the shift of to- (T) via the temperature dependence of 
D-term condensation. (The second effect was not con- 
sidered in Ref. 0.) In this paper, we have included both 
of these effects and calculated the density parameter of 
vr. 

In calculating the temperature-dependent mass of left- 
handed sneutrino m 2 L (T), one may think that it is also 
necessary to take account of the thermal mass of the sneu- 
trino. At the temperature T > Tc, the thermal mass is 
from the gauge-boson- loop diagrams, and is estimated as 



Sml L = l(3 ff 2 + 5 2 )T 2 . 



(40) 



Importantly, however, we are interested in the case where 
v(T) is non- vanishing since the production of vr is effec- 
tive in such a case. Then, the thermal mass of left-handed 
sneutrino should be suppressed since gauge bosons ac- 
quire masses from the expectation value of the Higgs bo- 
son. Detailed calculation of the thermal mass is more 
complicated in this case and we do not go into the de- 
tail. Instead, we perform our analysis with Eq. (|40[) and 
with Sm 2 L — (i.e., with neglecting the thermal mass 
of left-handed sneutrino). As we will see, production of 
right-handed sneutrino becomes effective at relatively low 
temperature and does not change much between two 
procedures. 





0.9 






0.8 






0.7 






0.6 




p 


0.5 




© 


0.4 






0.3 






0.2 
0.1 


1 AM 




-i i i i 1 1 1 1 1 1 i i i 1 1 1 1 i i i 1 1 i i i i rm 



20 40 



60 80 100 120 140 
T [GeV] 



FIG. 7: The left-right mixing angle of sneutrinos at finite 
temperatures 0(T) normalized bu O(0). We take here mc fl = 
100 GeV, and the mass difference AM = mp L — mp R is 1, 10 
and 100 GeV from left to right. The results obtained by using 
5ml, = and 5ml, given in Eq. (|40[) are shown by the solid 
and dashed lines, respectively. 



In Fig. [7] we show the left-right mixing angle at 
the finite temperature (relative to the zero-temperature 
value). We can see that the mixing angle is drastically 
suppressed when the temperature becomes close to Tc- 

Now, we estimate the relic density of vr coming from 
the decays of superparticle in the chemical equilibrium. 
We calculate J1~ E by solving Eq. (f30|) . In Fig. [3 we show 
il~ E /i 2 as a function of m;, . As we have mentioned, J1£ E 

Vr 1 Vr 

becomes larger as m VL gets close to rn UR . 7 In particu- 
lar, £1-^ becomes consistent with the density parameter 
of dark matter with the mass degeneracy of ~ 10 % for 
lay | < 3. We have also checked that fi? E decreases as 
increases. This is due to the fact that the sneutrino pro- 
duction at higher temperature becomes ineffective since 
the left-right mixing angle of sneutrinos is suppressed at 
high temperature. (See Fig. [9]) 

Eq. (|33|) suggests another possibility of enhancing £1 ~ E . 
By assuming degenerate neutrino mass, Yukawa coupling 
constants of neutrinos become larger and hence more vr 
can be produced. With Eq. (|33|) . for example, observed 
dark-matter density given in Eq. ([1]) is explained when 
neutrino masses are larger than 0.39 eV with degenerate 
neutrino masses. Notice that this value of the neutrino 
mass is consistent with the upper bound on the neutrino 
mass derived from the WMAP three-year data. 

Up to now, we have discussed the two sources of the 
Vr relic density, i.e., and J1? E . If both of them are 
sizable, the total present relic abundance is given by the 



7 Compared with the previous analysis in Ref. Q, Q9^h 2 is re- 
duced since the left-right mixing angle of sneutrinos becomes 
smaller due to the additional thermal effects. 
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FIG. 8: The relic density fig/i 2 of the LSP v R in terms of 
the left-handed sneutrino mass raa L . The results obtained by 
using 8m% L = and 5m% L given in Eq. (|40)l are shown by 
the solid and dashed lines, respectively. We take a„ = 1, 3, 5 
from left to right. Here ma R = 100 GeV, = 300 GeV and 
fiH = 150 GeV. The horizontal dot-dashed lines correspond 
to the dark-matter density ([l}. 
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FIG. 9: The relic density f2^/i 2 of the LSP v R in terms of 
the wino mass. The results obtained by using 5m% L = and 
Smp given in Eq. (|40[) are shown by the solid and dashed 
lines, respectively. We take here a v = 3, ma R = 100 GeV, 
mp L = 110 GeV and \ih = 150 GeV. The horizontal dot- 
dashed lines correspond to the dark-matter density {TJ. 

sum of these two contributions: 

f^ R =Qg+fig. (41) 

Before closing this subsection, we comment on the case 
where the right-handed neutrinos have Majorana mass. 
If the Majorana mass is smaller than the supersymme- 
try breaking masses of superparticles, vr may still be 
the LSP. In this case, however, the universe is likely to 



be overdosed with i>^-LSP. This is because y v becomes 
larger than the value given in Eq. with the Majorana 
mass. In this case, £lp E becomes much larger than the 
value we obtained for the purely Dirac case unless the 
reheating temperature is so low that superparticles can 
never be thermalized. (See also 0].) 



D. Other possibilities 

So far, we have discussed the production of vr by the 
decays of MSSM superparticles. In the rest of this sec- 
tion, we briefly comment on other possibilities. 

First, vr is potentially produced by the decays of non- 
MSSM superparticles. One interesting candidate is grav- 
itino. Since gravitino has only the gravitationally sup- 
pressed interactions, its lifetime is very long and gravitino 
of mass ~ 100 GeV decays after BBN epoch. When the 
primordial abundance of gravitino takes relevant value, 
its decay produces the correct amount of vr and the vr- 
CDM can be realized. In usual cases, decay of the grav- 
itino after BBN is severely constrained by the BBN con- 
straint 0,[2^|. However, if the gravitino is the NLSP and 
also if it decays only into vr + v (and v R + v), the BBN 
constraints are significantly relaxed. 8 This is because the 
branching ratios of gravitino decaying into hadrons and 
charged particles are very small in this case. We should 
note that, even in this case, the decay of the MSSM-LSP 
into gravitino is also constrained from BBN. To avoid it, 
the MSSM-LSP should dominantly decay into vr (and 
something else) with relatively small hadronic branching 
ratio or relatively short lifetime. 

In addition, in this paper, we have assumed that the 
initial abundance of vr after the reheating of inflation 
is negligible. Note, however, that we cannot exclude the 
possibility that vr is produced at the reheating (or also 
at the preheating) epoch by the decay of inflaton. This 
is possible when vr couples to the inflaton although it is 
highly dependent on the details of the inflation model. 

There is yet another possibility, vr might be produced 
as a coherent oscillation. This is because vr has an al- 
most flat potential and its initial amplitude after the in- 
flation might be displaced from the origin. In this case, 
if the initial amplitude is of the order of 10 9 GeV, such 
coherent mode of vr might be the dark-matter assuming 
that the reheating of the inflation completes before the 
start of the vr oscillation. 



8 Similarly, when gravitino is the LSP and ur is the NLSP, the 
gravitino dark matter may be possible without conflicting with 
the BBN observations. 
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IV. CONCLUSIONS 

We have considered the scenario where right-handed 
sneutrino becomes the CDM of the universe in the frame- 
work that neutrino masses are purely Dirac-type. We 
have investigated the production of vr by decays of su- 
perparticles in the early universe. Especially, the decay of 
the MSSM-LSP after its freeze-out time have been stud- 
ied in detail by using the minimal supergravity model of 
supersymmctry breaking. The NLSP is considered as the 
MSSM-LSP, xi or n. 

We have found that there is a wide parameter range in 
which vr is the LSP, and also that the relic density of the 
LSP vb coming from the late MSSM-LSP decay can be 
consistent with the dark-matter density in such a region. 
Assuming the minimal supergravity model, the mass of 
the vr is bounded not only from below but also from 
above. The lower bound comes from the phcnomcno- 
logical constraints (e.g., the Higgs boson mass and the 
rate of b — ► sj) or the requirement of vr being the LSP. 
On the other hand, the upper bound comes from the 
dark-matter density. In addition, in our framework, the 
MSSM-LSP is the MSSM-LSP is the Bino-like lightest 
neutralino (almost right-handed lighter stau) when m^ R 
is small (large). It should be noted that these features 
are insensitive to the neutrino mass hierarchy. 

We have also shown that the decay of the MSSM-LSP 
receives constraints from BBN, which leads to the upper 
bound on the mass of the MSSM-LSP. Indeed, it for- 
bids the MSSM-LSP with mass > 1 TeV in any case. 
It has been found that, when tan/3 = 10 and /in > 0, 
the X? MSSM-LSP can be consistent with the BBN even 
with the hierarchical neutrino masses, if the mass of Xx is 
small enough. On the other hand, the fx MSSM-LSP is 
severely constrained from BBN argument. We found that 
its decay can be consistent with the BBN only when the 
neutrino masses are degenerate and m„ is large. How- 



ever, the BBN bound strongly depends on the lifetime 
and becomes irrelevant when the lifetime is shorter than 
about 0.1 sec. If we go beyond the assumptions used in 
our analysis, such a short lifetime can be achieved. In 
such a case, the decay of the fx MSSM-LSP is cosmolog- 
ically harmless. 

When the CDM is indeed vr, it is not easy to test di- 
rectly its existence because of very small neutrino Yukawa 
coupling constants. However, the vr-CDM scenario has 
features quite its own, and we will have indirect hints 
from the future experiments and observations. Since 
neutrino masses are purely Dirac-type, experiments of 
the neutrinoless double-beta decay should give null re- 
sults. The direct searches for dark matter also should give 
null results. If supersymmetry would be found with the 
above-mentioned circumstance, vr would be regarded as 
a serious candidate of the LSP and CDM. In addition, the 
vr dark matter scenario may be distinguishable from the 
conventional scenario where the Xx MSSM-LSP is dark 
matter. One interesting feature of the vr-CDM is that 
the MSSM-LSP can be charged (or even colored), e.g., fx 
as we have considered. Even if the MSSM-LSP is Xii as 
we have shown, the mass relations between m^o and 
and also between and are significantly different 
in these two dark-matter scenarios. These differences will 
be testable in the precise measurements of superparticlcs 
at future collider experiments. 
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